Equilibrium temperature anisotropy and black-hole analogues 
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When long-range interactions are present the usual definition of temperature implies that two 
systems in thermal equilibrium can be at different temperatures. This local temperature has phys- 
ical significance, if the sub-systems cease to interact, each system will be at their different local 
temperatures. This is formally related to redshifting of temperature in general relativity. We pro- 
pose experiments to test this effect which are feasible using current microfabrication techniques. It 
is also possible to display thermodynamical analogues to black-hole space-time. 
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INTRODUCTION 

When a system possesses long-range interactions, 
many familiar notions of thermodynamics break down; 
the micro-canonical and canonical ensembles become in- 
equivalent ll[ , there may be no stable equilibrium config- 
uration [H, y], and heat-capacities can be negative [1,01 
as observed in fragmenting nuclei [5| and atomic clus- 
ters 0. Formally, standard thermodynamics has only 
been proved valid when interactions are short-range, or 
when the long-range interactions are screened (e.g. plas- 
mas) [7]. There are few methods for studying the ther- 
modynamics of systems with long-range interactions, al- 
though some models have been studied for special cases 
where the thermodynamic limit exists [8j], and under- 
standing them outside the standard framework using 
Tsallis or Renyi entropy has been attempted. 

Using a general formalism for studying such systems Q 
based on techniques used in general relativity (where the 
equivalence principle is exploited) - an effect was noted 
which we summarise and extend as follows. Consider 
many strongly interacting sub-systems in thermal equi- 
librium. Using the standard definition of temperature 
(defined as the local temperature), each sub-system is at 
a different temperature even though the entire system is 
at thermal equilibrium 0. Clearly the standard defini- 
tion does not satisfy a basic notion - that it be constant 
throughout the sample at equilibrium, yet it has a physi- 
cal meaning - if the interaction is turned off suddenly and 
the sub-systems isolated, they will be at their local tem- 
peratures. The observed temperature difference when a 
system is broken down into its parts is a property of the 
system and is a function of its self-interaction. What's 
more, the effect is formally related to effects found in 
curved space and black hole thermodynamics . Here we 
show that this anisotropy of local temperature can be ob- 
served, though this requires the long-range interaction to 
be strong compared with the characteristic temperature 
of the sample. However, recent advances in microfabrica- 
tion may allow experimental access to thermodynamical 
effects not found in macroscopic systems. An experiment 
to measure the break-down of temperature in quantum 



systems was recently proposed in EH- 

We review the basic formalism for analysing systems 
with long range interactions and then apply this to a pro- 
posed experiment. We also show that one can observe 
additional effects closely related to thermodynamical re- 
lations in black-hole space-times. 

LONG-RANGE INTERACTIONS 

Consider two sub-systems, with total energy m of the 
form 



m = E 1 +E 2 + G(E 1 ,E 2 ), 



(1) 



where Ei are the local (non-interacting) energies and 
G{E\,E 2 ) is some interaction potential (which may in- 
clude self-interacting terms). As an example, we consider 
the interaction of two clusters of classical spins in local 
magnetic fields b\ and b 2l with spin-spin interaction 

m = b\e\ + b 2 e 2 — J\e\/2 — J 2 e\/2 — J\ 2 e\e 2 , (2) 

where e\ and e 2 are the spin-excesses of each cluster = 
nj — n\, (nj and n\ being the number of up and down 
spins respectively of cluster i = 1,2), Ji are intra-cluster 
couplings, and J\ 2 is the inter-cluster coupling. Such an 
interaction arises from a standard interaction between 
spins of the form 



in 



E E 



J, 



(3) 



with the spins of cluster i being a\ . For small clusters, the 
interaction within each cluster can be approximately the 
same for all spins, not just nearest neighbour, i.e. J\\ = 
Ji > J 22 — J2 7 and the interaction between each cluster is 



approximately uniform over each cluster i.e. Jf^ 



We then drop constant terms to obtain Eq. ([2|). Gravity 
is another common example of such an interaction, where 
the potential term depends on the local energies (masses) 
of each system. 

One can define the temperature of one of the clusters 
(the system S), by making the other system very large 
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(the reservoir 1Z). If the total energy m is fixed, in this 
limit the probability that the system has energy Eg is 
given by [9( 

P (E s )dE s = ^s{E s )n R {E)e- Es ^dE s /Z m (4) 



where the sum is taken over all E = En + Eg consistent 
with total energy m = E + G(E$, En) and 



Z m — 



dE s n s (E s )e 



(5) 



m 



i.e. Z m is the total number of states at fixed total energy 
m. We then define the inverse temperature /3g in the 
usual manner in terms of the local extensive entropy 



fa 



dSn(E n ) 



dE- 



■R 



(6) 



We shall refer to (3e as the local temperature. The moti- 
vation for using this term comes from general relativity. 

Note that the temperature of the system is defined 
in terms of the derivative of the reservoir's entropy. In 
the non-interacting case, no issues arise from this defini- 
tion: if two systems are in thermal contact in the micro- 
canonical ensemble, then dS g E Es ^ — ^~§js ■ This is also 
true when the reservoir has no long-range interactions, 
or when the division of a single system into a reservoir 
and smaller system is purely formal (as we will see from 
symmetry considerations). In general however, this is not 
necessarily true - a point which will be discussed shortly. 
One therefore should keep in mind that the temperature 
is a property of a reservoir - it gives the distribution asso- 
ciated with a smaller system in contact with it. Finally, 
note that the local temperature, as defined, is a function 
of E, we make this explicit by writing /3g. There will 
be different "temperatures" depending on the value of E 
the entire system is found in. 

Let us now show that at equilibrium, a system can have 
an anisotropy in local temperature. Here we consider 
a slightly different situation and derivation than in 
For concreteness, one could consider a system, as above, 
composed of two clusters as above labelled 1 and 2. For 
fixed Ei and E 2 , the spin excesses of each cluster are 
fixed, thus their respective entropy is also fixed and given 
by 



Si(Ei) =-n\ log- 



,1 



,1 



nj log ■ 



,1 



(7) 



and since each entropy is determined only by the spin 
excess, the total entropy is additive, i.e. the number of 
states accessible to two systems is given by the number 
of states accessible to one system, times the number of 
states accessible to the other. I.e. 



n 12 (E 1 ,E 2 ) = a 1 (E 1 )n 2 (E 2 



(8) 



where O12 is the total number of states of the com- 
bined system, and Qi is the number of states accessi- 
ble to system i (the entropies of each system are just 
Si(Ei) = logfii(£?i)). Additivity of entropy holds be- 
cause we will find that we only need it at each value of 
Ei, E 2 and in this case, the entropy is additive in many 
systems because often the only correlations in a system 
are related to correlations in energy. We can thus con- 
sider more general systems with additive entropies (when 
conditioned on local energy) . We now consider these two 
systems in contact with a reservoir 1Z at fixed total en- 
ergy m = Ei + E 2 + E-jz + G(Ex,E 2 ) and again tak- 
ing the total local energy as E = Ei + E 2 + E R . The 
probability that the two systems have energy Ei and E 2 
for given local energy E is Qi(Ei)fl 2 (E 2 )$} R (E R )/Z E , 
and the probability that the local energy is E is de- 
noted by p{E) — ZE/Z mi with the partition functions 
z m = Y,e 1 e 2 e ^i{Ei)VL 2 (E 2 )VL R (E R ) at fixed m and 
z e = EEuE 2 '^i( E i)^(E 2 )n R (E R ) at fixed E. Then 
the probability that system 1 has energy Ei is 



Pi 



{E 1 )= Y,P(Ei,E 2 \E)Z E /Z„ 



E,E 2 



(9) 



ftl(£l)^2(£ 2 )ttfl(£fl)/Z m 

E,E 2 

Q±V!ll e SR{m)-p Ei ^ n 2 (E 2 )e : 



E,E 2 



where we have used the approximation that E-ji ;g> 
Ei + E 2 + G(Ei, E 2 ) to expand S R around m with the 
inverse global temperature defined as f3 = dS R (m) / dm. 
Since m is held fixed, an absence of interaction G{E\, E 2 ) 
would give the usual fact that pi (E± ) will be proportional 
to fli(Ei)e^ l3oEl , and symmetrically for p 2 (E 2 ) - thus 
the systems will be at equal temperature. Here, we have 
the additional sum over E and the factor due to the in- 
teraction term, which makes the distribution of system 1 
and system 2 asymmetric. Expanding G(Ei,E 2 ) around 
{E\,E 2 ) = (Ei,E 2 ) the average energies, gives 

0G{B 1 ,E 2 ) , 



p(£?i) oc ni{Ei)e 



-(n 



(10) 



which looks like a canonical distribution. By symmetry, 
each system thus behaves as if it has an inverse temper- 
ature 



(ii) 



One can calculate that this local temperature matches 
the standard definition of temperature given by Eq. © 
and we thus see that there is a temperature anisotropy 
whenever dG{E 1 ,E 2 )/8E 2 ± dG(E u E 2 )/dE 1 

Let us now understand the physical meaning of this 
local temperature. Eq. {4| gives the probability dis- 
tribution of a system S in terms of the local temper- 
ature /?, and the local energy Eg. However, the lo- 
cal energy is not a conserved quantity, and does not 



3 



contain the interacting term. To see what the locally 
conserved energy is, we can expand m to first order 
m(E 1 ,E 2 )~m(E 1 ,E 2 ) + §§ [ (E 1 -E 1 ) + §£(E 2 -E 2 ). 

Since constant terms can be ignored, E l a = §^-Ei can be 
identified with the energy of each system in the presence 
of the mean field due to the other interacting systems 
and is (1 + aG ^'" Ea ^ )Ei. In the case of a single system 

and reservoir, this energy is E = Es-§^ and serves to 
tell us how we should define the locally conserved energy. 
As is expected for a non-extensive system, the quantity 
is also non-additive, but, as with the local temperature, 
will correspond with what happens in general relativity. 

E a is also a function of E but we will not write this 
explicitly. With respect to energy levels E Ql each sys- 
tem acts approximately as if it is at inverse temperature 
P (although in actual fact it is a mixture of tempera- 
tures @). We can think of E Q as the effective energy, 
i.e. it is the energy of the system in the presence of in- 
teraction with another system, and j3 Q can be thought of 
as the closest thing one has to a physical temperature - 
looking at Eq. (fTTj) we see that the global inverse temper- 
ature (3 will be equal for both systems. Note however, 
that the global temperature is not an intensive quantity. 

If the interaction term is suddenly "turned off", (one 
can imagine that the spins are suddenly separated so 
that they no longer interact), then the energy that is 
ascribed to each energy level is no longer E l , but rather 
its local energy E\ or E 2 . If the change happens ex- 
tremely quickly then the overall state of the system will 
not change. Since 

fcEi = p E l (12) 

and the new energy of each system is now Ei, the mea- 
sured temperatures will be /Jj - the local temperature. In 
these systems, the local temperature becomes the physi- 
cal one when the interaction is turned off, which is anal- 
ogous to the fact that in general relativity, the local tem- 
perature is the physical one measured by free-falling ob- 
servers (for whom the gravitational interaction is "turned 
off"). 

Using the two coupled Ising models of Eq. © we would 
get a temperature difference of 

1 - JuPi/h - Jizilb 2 



(h = Pi 



(13) 



1 - J\ 2 e 2 /bi - .hei/bi 

Now if initially these two systems (or clusters), are 
far apart, and at equal global temperature, one cannot 
push them together both adiabatically and isothermally 
(constant global temperature) as is possible in the non- 
interacting case. This can be seen from Eq. (fT2"|) . Moving 
the systems together adiabatically requires keeping E (3 
fixed, but since E a changes when J\ 2 becomes significant, 
one cannot keep [3 constant. By recalculating E D one can 
therefore calculate the new global temperature. We see 
therefore that the global temperature is not an intensive 
quantity. 



BLACK-HOLE ON A BENCHTOP 

To obtain an analogue of a black-hole, we consider a 
single cluster of spins with long-range self-interaction. 
This can be obtained by setting h 2 = J 2 = J\ 2 = 0. This 
gives a relation between the local and global temperature 
of 



P=±Po\ l 



2J 2 m 



(14) 



The positive solution looks exactly like the relation- 
ship between global and local temperature in the 
Schwarzschild black-hole space-time due to redshifting, 
if we set the radius r = h 2 and equate the spin-spin 
coupling with Newton's constant Ji = G. Here (3 a is 
the temperature at infinity, while f3 would be the local 
temperature measured by an observer sitting close to the 
black-hole horizon. In the analogue, the horizon is real - 
for fixed m, there is no value of spin-excess which allows 
r < 2Gm i.e. m > h 2 /2J\. There are also many similar- 
ities between this analogue and a black hole in terms of 
the way the entropy and energy behave, and we refer the 
reader to d, [ijj for details. 



EXPERIMENTAL REALISATION 

Ideally, we would like to study these phenomena in 
complex macroscopic systems. However these usually 
do not possess a strong enough long-range interaction 
to produce an appreciable effect, hence it is necessary to 
go to small systems where the the coupling between sub- 
systems can be comparable to the local energies and low 
temperatures are easier to achieve. 

The simplest system to study that of two spin clusters 
coupled via an Ising-like interaction, 



H = bi ■ a 1 



- #2 i 1 



(15) 



where bi are local magnetic fields and J is the interac- 
tion energy between the two clusters. We can map this 
system to a pair of pseudo-spins such as double quantum 
dots with a excess electrons whose localisation in one or 
the other dot defines the pseudo-spin vector (Fig. [1] or 
Fig. [5]). A single electron transistor measures the total 
charge excess in one dot or the other. 

To implement the separation of the two sub-systems 
before measurement, the device in Fig. [3] could be used. 
A surface acoustic wave (SAW), produced in a piezoelec- 
tric material by an inter-digital transducer driven by a 
microwave generator, induces a travelling sinusoidal elec- 
tric field in a 2-dimensional electron gas (2DEG) situ- 
ated just below the surface of a modulation-doped GaAs- 
AlGaAs semiconductor. Bias electrodes on the surface of 
the device deplete the 2DEG of conduction electrons on 
the region of the system. 
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FIG. 1: Coupled double quantum dot (DQD) system con- 
structed from donor phosphorus atoms in silicon with 20nm 
vertical and 30nra lateral separation. Each left and right 
DQD pair is initialised with one electron shared between the 
two dots. The local Hamiltonians can be controlled by elec- 
trodes affecting the tunnelling rate and bias between the two 
dots. The interaction is of the Ising type (J = 2.9meV) and 
single electron transistors (SETs) measure the localisation of 
the electrons to the upper or lower dots of each pair. 
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FIG. 2: Surface acoustic wave (SAW) device. The SAW cre- 
ated by an interdigital transducer produces a sinusoidal piezo- 
electric field in the direction of propagation which, together 
with surface electrode located above a 2DEG, create a pattern 
of moving quantum wells which can be manipulated by static 
voltages. The advantage of this configuration is that the in- 
teracting sub-systems can be separated before measurement 
of their respective populations. Correlated current readout 
measures the population in the position basis. 



By suitable electrode potentials at the entrance of the 
system, an exact number of electrons can be transported 
in moving quantum dots, along quasi-lD quantum chan- 
nels (denned by surface electrodes) through the device. 
Double well potentials can be created with a specified 
number of electrons in each, and by using suitable elec- 
trode geometry, made to interact via an Ising coupling. 
Travelling along a sufficiently long channel, the sub- 
systems arc allowed to thcrmalise, after which the two 
sub-systems are rapidly separated and read-out by sim- 
ply measuring the current via Ohmic contacts. 
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FIG. 3: A plot of the temperatures of two interacting 
sub-systems consisting of several electrons with Coulomb ex- 
change energy. 
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